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Abstract 

We investigate the behavior of single polyelectrolytes in multivalent salt solutions under the 
action of electric fields through computer simulations. The chain is unfolded in a strong electric 
field and aligned parallel to the field direction, and the chain size shows a sigmoidal transition. 
The unfolding electric field E* depends on the salt concentration and scales as V^ 1 ^ 2 with V being 
the ellipsoidal volume occupied by the chain. The magnitude of the electrophoretic mobility of 
chain drastically increases during the unfolding. The fact that E* depends on the chain length 
provides a plausible mechanism to separate long charged homopolymers by size in free solution 
electrophoresis via the unfolding transition of globule polyelectrolytes condensed by multivalent 
salt. 

PACS numbers: 82.35.Rs, 87.15.Tt, 36.20.Ey, 87.15.ap 
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Electrophoresis is an important separation technique in molecular biology and chemical 
analysis [l, 2]. Nonetheless, in a free solution, polyelectrolytes (PEs) like DNA molecules 
cannot be size-separated in an electric field [3j owing to the proportionality of the hydro- 
dynamic friction and the total charge of the molecule to its chain length l|, |4| . Therefore, 
to separate homo-PEs by size is commonly carried out in a sieving matrix, like gel |2|, in 
which biased reptation is responsible of the separation mechanism . However, efforts 
have been continually put into free solution electrophoresis (FSE) due to its high throughput 
and analysis efficiency, compared to the gel electrophoresis (GE). A breakthrough was done 
by a method called "end-labeled free-solution electrophoresis" (ELFSE) 6]. In ELFSE, the 
friction to charge balance is broken by attaching uncharged giant molecules orpolymers to 
the ends of the PEs, producing a length-dependent electrophoretic mobility 3). Recently, 
Netz proposed a new strategy to separate homo-PEs by size based upon drastic change of 
migration speed in the unfolding of collapsed PEs by an electric field [8!]. He predicted 
that the electric field to unfold a chain, E* , follows the scaling law N~ 3u ^ 2 where N is the 
chain length and v the swelling exponent. Therefore, longer chains will be unfolded and 
separated out using a smaller electric field. It could be a plausible way to sequence DNA, 
providing a read- length of hundred-kilo bases [8j, in comparison with thousand bases in 
capillary GE and hundred bases in ELFSE [7|. In this study, we use a coarse-grained model 
to verify the prediction of Netz for the first time by means of computer simulations. The 
objective is to understand how electric fields influence the conformation of single chains and 
the electrophoretic mobility. 

There are many ways to collapse PE chains 9j, for e.g., reducing the temperature, de- 
creasing the dielectric constant of solvent, modifying the chain charge density, and so on. 
Adding multivalent salt into the solution is a simple way to initiate the collapse at ambient 
condition. It is the condensation of multivalent counterions on the chains which induces 
the chain collapse [10]. The collapsed PEs exhibit the most compact structure at the salt 
concentration C* where the multivalent counterions and the chains are in charge equiva- 
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111 ] , re sulting in a nearly zero effective chain charge and hence zero electrophoretic 
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131 ] . The condensation of multivalent counterions also leads to two very relevant 



phenomena, overcharging and charge inversion 12j. We showed in a recent study [13| that 
although charge overcompensation (i.e. overcharging) occurs on the surface of PE at high 
salt concentrations, the effective chain charge may not invert its sign; in other words, charge 
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inversion does not necessarily happen with overcharging. In that study, the effective chain 
charge was calculated by the ratio of the electrophoretic mobility to the diffusive mobility, 
and the applied electric field was very weak, keeping unmodified the chain conformation. 
The fact that a strong electric field could largely change the conformation and mobility 
lightens up a possible way to separate PEs by size and drives our motivation to study the 
electrophoretic behavior of PEs in multivalent salt solutions. This kind of study, interplay- 
ing at the same time the strength of electric field and the salt concentration, is still few 
in literatures, specially by a method able to provide information at molecular level. It will 
enrich our understanding of the behavior of PEs in the presence of multivalent counterions 

V 

in a whole range of electric field. In Netz's study [8J, a PE was collapsed in a salt-free 
solution by decreasing the dielectric constant. Although his results gave a valuable picture 
of PE behavior, a more realistic system such as salt-collapsed PEs in electric fields should 
be investigated in detail in order to put into reality of his idea. Moreover, to unravel a 
DNA molecule is a key issue in DNA sequencing and diagnostics. It can be done by several 
methods such as optical tweezers [14j], flows [15], and electric fields 16[. This work can 
provide deep insight of the unraveling of charged biopolymers in electric fields. 

We simulate a single PE chain using a bead-spring chain model, consisting of N 
monomers. Each monomer carries a negative unit charge — e and the adjacent monomers 
are connected by virtual springs described by a finitely extensible nonlinear elastic poten- 
tial U sp (b) = — 0.5fc6^ aa . ln(l — (& 2 /&max)) where b is the bond length, b max the maximum 
extention, and k the spring constant. We add (4:l)-salt into the system and hence there 
are three kinds of ions in the solution: the tetravalent cations (counterions) and monovalent 
anions (coions) dissociated from the added salt, and the N monovalent cations dissociated 
from the chain. The excluded volume of the monomers and ions is modeled by a shifted 
Lennard- Jones potential U^j(r) = e^j [2(a/r) e — l] 2 truncated at its minimum where a rep- 
resents the particle diameter and £lj the coupling strength. The particles also interact 
with each other via Coulomb interaction U e i(r) = kBTXsZiZj/r where ks is the Boltz- 
mann constant, T the absolute temperature, Z± the charge valence, r the distance and 
Xb = e 2 / {ATieQeksT) the Bjerrum length (eo the vacuum permittivity). The solvent is con- 
sidered as a continuum medium of dielectric constant e and its effect on the chain and ions 
is implicitly taken into account by Langevin equation irtifi = —Qfi + F c + ZieEx + fji where 
rrii is the ith particle mass and Q the friction coefficient, F c = —dll/dfi the conservative 
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force acting on the particle, and ffi the random force satisfying the fluctuation-dissipation 
theorem: (f)i(t) ■ f)j(tf)) = 6kBTQ5ij5(t — t'). The system is placed in a periodic box, sub- 
ject to an external electric field Ex pointing toward x-direction. Ewald sum is applied to 
calculate Coulomb interaction. We know that in dilute solutions without electric fields, the 
dynamics of chain is described by Zimm dynamics [171 ] . However, under the action of electric 
fields, the coions and the counterions move in opposite directions, which largely cancels out 
the hydrodynamic effect in the solutions. As well documented in the literatures 0, [l8 |. 
hydrodynamic interaction is shielded in a typical electrophoretic condition. Therefore, we 
neglect it in this study. We remark that this effect will become important when the chain 
length is very short 19|, |20j. 

We assume that all the monomers and ions have identical mass m and diameter a. In the 
following text, we choose a, m, ksT to be the unit of length, mass, and energy, respectively. 



Hence, the time unit will be crym/ '(/cgT), the unit of concentration a~ 3 , the unit of electric 
field ksT/^ea), and so on. We set Elj = 0.8333, k = 5.8333, b max = 2, A# = 3, Q = 1. 
The chain length N is varied from 12 to 384 and the volume of simulation box is changed 
accordingly to keep constant the monomer concentration C m at 3 x 1CT 4 . For a short chain 
(N < 48), we use a cubic box whereas for a long chain (N > 48), a rectangular box elongated 
in the field direction is used to prevent chain overlapping under periodic boundary condition. 
The sides of the rectangular box are equal to 54.288 in y- and z-directions and 54.288iV/48 
in x-direction. Salt concentration C s is varied from 0.0 to 6 x 10 -4 . The strength of the 
electric field E is changed from 0.0 to 2.0. We launch Langevin dynamics simulations 21] 
with a time step equal to 0.005. A pre-run of 5 x 10 6 time steps is performed to bring the 
system to a steady state and a production run is then executed for, at least, 5 x 10 7 time 
steps. 

In Fig. [T^a), we show the variation of chain size for a PE of N = 48, characterized by the 
square of radius of gyration, R 2 g1 against C s under the action of different electric field E. In 
zero field, the chain size decreases firstly, due to the screening or the bridging of tetravalent 
salt, up to C s = C* = C m /A{= 7.5 x 10~ 5 ), and then progressively increases. This size 
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22 1 . Because we focus here on the 



decrease is the well-known salt-induced condensation 
behavior of condensed PEs in electric fields, the value of C s is not so elevated to be able 
to observe the chain decondensation. At a really high C s , the PE is expected to attain the 
size of its neutral counterpart 23| . Under the action of an electric field, R g becomes larger. 
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FIG. 1: (a) R 2 g vs. C s in different E. (b) R e /L c vs. E at different values (xl0~ 5 ) of C s . The "x" 
symbols denote the E* at a given C s plotted on the associated curve. The error of data in this 
Letter is either smaller than the symbol size of data or indicated directly by the error bar. 



The stronger the electric field, the larger the deviation from the chain size in the zero field. 
Moreover, there is a threshold of electric field near E = 0.4, smaller than which the chain is 
basically not deformed by the field at C s > C*. It shows that a coil PE (in C s < C*) is easier 
to be deformed than a condensed PE (in C s > C*). For E > 0.4, the chain reexpansion is 
suppressed by the strong field and R g decreases with C s . A very strong field (E = 2.0) can 
completely prevent the collapse of the PE and R g is a constant, indicating a chain unfolding. 
We have verified that the unfolded chain is aligned parallel to the field direction. The degree 
of unfolding can be studied by calculating the ratio of the end-to-end distance R e to the 
contour length L c = (N — 1)6 and the results are plotted in Fig. [Tfb). Each curve denotes a 
salt concentration and exhibits a sigmoidal variation against E. The sharp transition of the 
curve suggests the existence of a critical electric field E* to unfold a chain, over which the 
chain is stretched and R e /L c drastically increases. In this study, the chain is unfolded up 
to a value of R e , 80% of the contour length, and hence displays like a rod. The transition 
occurs at different E, showing that E* depends on C s . 
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We determine E* by equating the polarization energy W po i = p ■ E/2 to the thermal 
energy /c#T where p is the induced dipole moment in the electric field E [8(. It is because 
the thermal fluctuation can destroy the polarization order while W po i < k^T and hence, E 
is too weak to perturb the chain orientation and conformation. In solutions, charged chains 
can form complexes with ions. To calculate p, we need to know first the "members" of the 
PE complex. We define the members to be the monomers and ions lying in a worm-shaped 
tube around PE, which is the union of spheres of radius r t , centered at each center of the 
monomers. We choose r t = 2. Inside this tube, the kinetic energy (3fceT/2) of a counterion 
is smaller than the Coulomb attraction to a monomer and hence, the ion is "bound" to the 
chain. We calculate p by the equation p = J2i Zie(fl — r cm ) where f\ is the position vector 
running over all the members of the complex and f cm is the center of mass of the PE. Fig. [2(a) 
shows the x-component p x of p as a function of E. We find that the linear response theory, 
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FIG. 2: (a) P x vs. E at different C s . The dashed line denotes the equation p x = 2ksT/E. (b) E* 
as a function of C«. 



p x = aE, holds true when E is weak, and breaks down soon after E surpasses 2ksT '/p x 
(plotted in dashed line). E* is determined by the intersection of the curve with the dashed 
line and is a function of C s , shown in Fig. [2(b). E* increases with C s , attains a maximum 
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value, and then progressively decreases. It shows that the most compact structure occurs 
at C s = C*. By plotting E* on the corresponding curve in Fig. QJb) (with the "x" symbol), 
we show that the E* denotes the onset where the PE starts to unfold. 

Netz has predicted E* ~ J\f~ 3u / 2 [sj but this scaling law has not yet been verified. Here, 
we study the dependence of E* on chain length by varying N from 12 to 384. We focus on 
two cases: C s = 0.0 and C s = C*. The results are presented in Fig. [3]^a). In the log-log plot, 




FIG. 3: (a) E* vs. N, (b) R 2 vs. N, and (c) (R^R^R 2 ^) 1 / 2 vs. N, at C s = 0.0 and C*. 



E* lies on a straight line, showing a power-law dependence. Effect of finite chain length is 
not important. Least square fits of the data (omitting the N = 12 point) yield jV" -1 - 05 ® 
for C s = 0.0 and jV -0 - 453 ^ for C s = C*. According to the Netz's theory, the swelling 
exponent should be 0.70(2) and 0.302(2) for the two cases, respectively. However, a direct 
study of the scaling law of R 2 in the zero field (Fig. [3](b)) shows v = 0.91(1) for C s = 0.0, 
of which Netz's prediction has large deviation, and 0.334(1) for C s = C*, with which it 
agrees within 10% error. We rederive here the scaling law of E*. The general polarizability 
is a tensor and can be written as a = €on~ l V where n is the depolarization tensor and 
V is the volume of the dielectric object [24] . For a fixed geometry, thus a fixed n, the 
polarizability is proportional to V. In this study, k B T = p x E*/2 = aE* 2 /2, which implies 

7 



E* ~ a -1 / 2 ~ V~ x l 2 . In Netz's theory, V was simply estimated by R^. Here, we do it in a 
more precise way, from the gyration tensor of the chain. The three eigenvalues of the tensor 
denote the squares of the gyration radii, R 2 gl1 R g2 , and R 2 3 , for the three principal axes of 
inertia. By considering the PE complex as an ellipsoid, we estimate V as {R 2 g iR 2 g 2R%) 1 ^ '• 
Fig. He) shows that {R 2 gl R 2 g2 R%) 1/2 follows a priori a power law. The effect of finite chain 
length is more important in this case because an ellipsoidal description is good only when 
chain length is long. Fitting the data for iV > 96 gives V ~ TV 216 ® and N im ^ for the 
two cases, which yields E* ~ N' 1 - 08 ^ and N" - 50 ^, in good agreement with the results of 
Fig. EKa). 

Knowing that E* depends on N, we study now the electrophoretic mobility of PE, fi pe , 
for different chain length. fi pe is calculated by v pe /E where v pe is the velocity of the center 
of mass of PE in the field direction. We focus on the case C s = C* where the chain exhibits 
the most compact structure. Fig. Hta) and (b) show, respectively, how /i pe and R e /L c vary 
with E. In weak fields, /i pe is zero, showing that the PEs are of charge neutral in the salt 




FIG. 4: (a) /j, pe vs. E and (b) R e /L c vs. E, for different chain length at C s = C*. (c) number of 
condensed tetravalent counterions on each monomer of the PE of N = 384 in different strength of 
E. 
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solutions. This confirms the experimental observations 10J. As soon as E is strong enough 
to unfold the chain, fi pe shows a drastic decrease in accompany with chain unfolding. The 
net chain charge is hence negative. The longer the chain length, the sharper the unfolding 
transition and also the decrease of /i pe . The fact that the occurrence of the // pe -decrease is N- 
dependent provides a good mechanism to separate PEs by the chain length. One advantage 
to do electrophoresis at C* is that only the chains being unfolded will be moved out by the 
electric field, leaving the shorter chains in the center. For N > 96, fi pe shows a plateau 
when the chain is completely unfolded, which is basically independent of chain length. This 
predicts a field-independent free-draining behavior when the chains are stretched. We have 
verified that in this region, the number of the condensed tetravalent counterions is roughly 
a constant. Fig. H^c) shows the variation of these ions distributed along the chain backbone 
of N — 384 in different E fields. The distribution changes from a horizontal line (of value 
0.25) at E = 0.18 where the chain is in a collapsed state, to a line inclined toward the field 
direction at E = 0.2 where the chain is unfolded, followed by a tangential (E = 0.4) and 
then a horizontal (E = 1.0) sigmoidal curve. In a high field region E > 1, the baseline of 
the distribution moves downward, showing a decrease of the number of the condensed ions, 
which implies the further decrease of fi pe . We have calculated the electrophoretic mobility 
of the condensed tetravalent counterions, /z + 4. The value is zero when the chain is collapsed. 
Once the chain is unfolded, /i + 4 becomes positive and increases with E; thus, the ions are 
not tightly bound to the chain in the electric fields but glide along the chain following the 
field direction. They can fly away off the chain and other ions condense onto it, establishing 
a dynamic equilibrium which keeps constant the number of the condensed ions. 

The results clearly demonstrate the possibility to separate PEs by size in multivalent 
salt solutions under the action of appropriate electric fields. Mapping our models to a PE 
in an aqueous solution at room temperature yields a = 2.38A and hence, the critical field 
E* = 0.28 for N = 48 is equivalent to 2.9 x 10 5 V/cm in real unit. Therefore, to unfold a 
PE chain with N = 10 6 demands E* of a value roughly 2 kV/cm according to the scaling 
law, which is feasible in experiments and can be applied to separate long polymer chains 
in compliment to the current FSE methods. Our results confirm the idea of Netz [8]. For 
practical reason, we wish the value of E* to be as small as possible. To achieve this goal, 
we can use trivalent salt, or even divalent salt, to condense PEs to a less-compact collapsed 
chain structure [Hj, which will reduce E*. We remark that multivalent salt can also induce 
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multi-chain aggregation. To diminish this effect, salt of small ion size should be used as the 
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23] 



condensing agent, which induces mainly single chain collapse 

In summary, we have studied the static and dynamic properties of homo-PEs in multi- 
valent salt solutions under the action of electric fields. The electrophoretic mobility of PE 
shows a drastic decrease due to the change of the effective chain charge while the chain 
is unfolded. The unfolding electric field E* sensitively depends on the chain length. By 
studying the scaling law of E*, we presented the first, strong evidence of the feasibility to 
separate homo-PEs by size, via unfolding transition, in free solution electrophoresis. 

This material is based upon work supported by the National Science Council, the Repub- 
lic of China, under Grant No. NSC 94-2112-M-007-023 and NSC 95-2112-M-007-025-MY2. 
Computing resources are supported by the National Center for High-performance Comput- 
ing. 
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counterion /i + 4 and that of chain function of E, is available free of charge via the 



Internet at http://pubs.acs.org 



* E-mail: pyhsiao@ess.nthu.edu.tw 
[1] J.-L. Viovy, Rev. Mod. Phys. 72, 813 (2000). 
[2] K. Kleparmk and P. Bocek, Chem. Rev. 107, 5279 (2007). 
[3] B. M. Olivera, P. Baine, and N. Davidson, Biopolymers 2, 245 (1964). 
[4] D. Long, J.-L. Viovy, and A. Ajdari, Phys. Rev. Lett. 76, 3858 (1996). 
[5] H. Cottet, P. Gareil, and J.-L. Viovy, Electrophoresis 19, 2151 (1998). 
[6] P. Mayer, G. W. Slater, and G. Drouin, Anal. Chem. 66, 1777 (1994). 

[7] G. Heller, G. W. Slater, P. Mayer, N. Dovichi, D. Pinto, J.-L. Viovy, and G. Drouin, J. 
Chromatogr. A 806, 113 (1998); H. Ren, A. E. Karger, F. Oaks, S. Menchen, G. W. Slater, 
and G. Drouin, Electrophoresis 20, 2501 (1999); R. J. Meagher, J.-L Won, L. C. McCormick, 
S. Nedelcu, M. M. Bertrand, J. L. Bertram, G. Drouin, A. E. Barron, and G. W. Slater, 



10 



Electrophoresis 26, 331 (2005). 

R. R. Netz, Phys. Rev. Lett. 90, 128104 (2003); J. Phys. Chem. B 107, 8208 (2003). 

S. K. Tripathy, J. Kumar, and H. S. Nalwa, eds., Handbook of Poly electrolytes and Their 

Applications, vol. I, II, III (American Scientific, Stevenson Ranch, CA, 2002). 

V. A. Bloomfield, Curr. Opin. Struct. Biol. 6, 334 (1996). 

P.-Y. Hsiao and E. Luijten, Phys. Rev. Lett. 97, 148301 (2006). 

A. Y. Grosberg, T. T. Nguyen, and B. I. Shklovskii, Rev. Mod. Phys. 74, 329 (2002); 
M. Quesada-Perez, E. Gonzalez- Tovar, A. Martin-Molina, M. Lozada-Cassou, and R. Hidalgo- 
Alvarez, Chem. Phys. Chem. 4, 234 (2003). 
P.-Y. Hsiao, J. Phys. Chem. B 112, 7347 (2008). 
S. B. Smith, Y. Cui, and C. Bustamante, Science 271, 795 (1996). 
T. T. Perkins, D. E. Smith, R. G. Larson, and S. Chu, Science 268, 83 (1995). 
S. Ferree and H. W. Blanch, Biophys. J. 85, 2539 (2003). 

M. Doi and S. F. Edwards, The theory of polymer dynamics, (Oxford, New York, 1987). 

M. Tanaka and A. Y. Grosberg, Eur. Phys. J. E 7, 371 (2002); M. Tanaka, Phys. Rev. E 68, 

061501 (2003). 

E. Stellwagen, Y. Lu, and N. C. Stellwagen, Biochem. 42, 11745 (2003). 

K. Grass, U. Bohme, U. Scheler, H. Cottet, and C. Holm, Phys. Rev. Lett. 100, 096104 (2008). 



The simulations were run using modified LAMMPS package (http:/ /lammps.sandia.gov/). 
M. Olvera de la Cruz, L. Belloni, M. Delsanti, J. P. Dalbiez, O. Spalla, and M. Drifford, J. 
Chem. Phys. 103, 5781 (1995). 

P.-Y. Hsiao, J. Chem. Phys. 124, 044904 (2006); Macromolecules 20, 7125 (2006). 

L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrodynamics of Continuous Media 

(Pergamon, New York, 1984), 2nd ed. 



11 



